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Abstract 

The singlet scalar model is a minimal extension of the Standard Model 
that can explain the dark matter. We point out that in this model the 
dark matter constraint can be satisfied not only in the already considered 
WIMP regime but also, for much smaller couplings, in the Feebly Inter- 
acting Massive Particle (FIMP) regime. In it, dark matter particles are 
slowly produced in the early Universe but are never abundant enough to 
reach thermal equilibrium or annihilate among themselves. This alter- 
native framework is as simple and predictive as the WIMP scenario but 
it gives rise to a completely different dark matter phenomenology. After 
reviewing the calculation of the dark matter relic density in the FIMP 
regime, we study in detail the evolution of the dark matter abundance 
in the early Universe and the predicted relic density as a function of the 
parameters of the model. A new dark matter compatible region of the 
singlet model is identified, featuring couplings of order 10" 11 to 1CP 12 for 
singlet masses in the GeV to TeV range. As a consequence, no signals at 
direct or indirect detection experiments are expected. The relevance of 
this new viable region for the correct interpretation of recent experimental 
bounds is emphasized. 



1 Introduction 

Weakly Interacting Massive Particles (WIMPs) constitute the most studied class 
of dark matter candidates. They appear in several extensions of the Standard 
Model and have the virtue of acquiring, via thermal freeze out in the early 
Universe, the right relic density to explain the dark matter. Among them, we 
can mention the neutralino in supersymmetric models PQ, the lightest Kaluza- 
Klein particle in Universal Extra-dimensional models [2J, as well as the singlet 
[3J S] and the doublet [5J IE1 [Z] scalar extensions of the Standard Model. At 
present, in spite of its popularity, the WIMP framework for dark matter has 
to be regarded as an interesting hypothesis to be tested in current and future 
experiments. WIMPs typically give rise to a number of signatures in accelerator 
searches and in direct and indirect dark matter detection experiments. In fact, 
it has been recently claimed that the moment of truth for WIMPs has finally 
arrived [5]; if they exist they will be discovered in the next 5 to 10 years. So 
far, however, there is no experimental evidence that dark matter is actually 
composed of WIMPs. 
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Moreover, viable and well-motivated alternatives to the WIMP framework do 
exist [HI EH E] ■ In [TT] , a new calculable mechanism of dark matter production, 
dubbed freeze in, was introduced. It involves a Feebly Interacting Massive 
Particle (FIMP) which never attains thermal equilibrium in the early Universe. 
FIMPs are slowly produced by collisions or decays of particles in the thermal 
plasma and, in contrast to WIMPs, are never abundant enough to annihilate 
among themselves. Thus, the abundance of FIMPs freezes in with a yield that 
increases with the interaction strength between them and the thermal bath. 
FIMPs, in addition, give rise to completely different signatures in colliders and 
dark matter experiments [11J . Due to its simplicity and predictivity, the freeze 
in of FIMPs is a strong contender to the WIMP framework for dark matter. 

The singlet scalar model [3J |2] is a minimal extension of the Standard Model 
that can account for the dark matter. It contains just one additional field, an 
scalar that is singlet under the Standard Model gauge group and odd under a 
new Zi symmetry, and two new parameters. The phenomenology of this model 
has been extensively studied in a number of previous works (see e.g. [12l IT3] ) . 
All of them, however, have assumed the WIMP solution to the dark matter 
constraint. In this paper, we point out that the singlet model can also satisfy 
the relic density constraint in the FIMP regime, and we analyze in detail the 
new viable region corresponding to that solution. After briefly introducing the 
model, we review the calculation of the relic density in the FIMP regime of the 
singlet model. Then, we study the evolution of the dark matter abundance and 
the predicted relic density as a function of the parameters of the model. The 
new viable region corresponding to FIMP dark matter is identified -it features 
couplings between 10 -11 and 10 -12 for singlet masses in the GeV to TeV range. 
Finally, we discuss the implications of these results for the correct interpretation 
of recent experimental bounds. 

2 The singlet scalar model of dark matter 

The singlet scalar model [3J 0] is one of the minimal extensions of the Standard 
Model that can explain the dark matter. It contains an additional field, S, that 
is singlet under the SM gauge group and odd under a new Z2 symmetry that 
guarantees its stability. The Lagrangian that describes this model is 

C = C SM + \d»Sd»S - ^S 2 - ?fs 4 - XS 2 H^H , (1) 

where Csm denotes the Standard Model Lagrangian and H is the higgs doublet. 
The above is the most general renormalizable Lagrangian that is compatible 
with the assumed symmetries and field content of the model. The singlet scalar 
model introduces, therefore, only two new parameter^: the singlet mass (7715 = 
y/rriQ + Xv 2 ) and the coupling to the higgs boson, A. In addition, the higgs 
mass, a SM parameter, also affects the phenomenology of the model. 

In the early Universe, singlets are pair-produced as the particles in the ther- 
mal plasma scatter off each other. The dominant production processes are s- 
channel higgs boson mediated diagrams originating in a variety of initial states: 

1 \g is essentially irrelevant. 
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//, W + W~, Z°Z°, and M@. Likewise, they can also be produced from the 
initial state hh either directly or through singlet exchange. As a general rule, it 
is the initial state W + W~ that tends to dominate the total production rate of 
dark matter in this model. 

The main advantage of the singlet model with respect to other models of 
dark matter is its simplicity, which allows one to make concrete predictions 
about dark matter observables. In fact, by imposing the relic density constraint, 
fls = Qdm, one can eliminate A and compute the expected rates for direct and 
indirect detection experiments as a function of the singlet mass only (for a given 
higgs mass). This program has indeed been implemented in a number of recent 
works -see e.g. [T2l [13]. In them, it was demonstrated that the singlet model 
can be probed in current direct and indirect detection experiments. 

All previous works on the singlet scalar model have assumed the WIMP solu- 
tion to the dark matter constraint. It has recently been realized, however, that 
an alternative mechanism of dark matter production, the freeze in of FIMPs, 
which is as simple and predictive as the freeze out of WIMPs, can also explain 
naturally the observed value of the dark matter density. The singlet scalar 
model has the particularity that both solutions to the dark matter constraint, 
WIMP and FIMP, can be realized -although in very different regions of the 
parameter space. In the following, we will focus on the FIMP regime of the 
singlet scalar model. 

3 The relic density in the FIMP regime 

It is easy to convince oneself that not only in the singlet model but in all 
models where the interactions of the dark matter particle are determined by a 
free parameter A, the dark matter constraint can be satisfied not only in the 
usual WIMP regime, where A typically lies between 1CP 1 and 10~ 3 , but also for 
much smaller couplings -in the FIMP regime. Imagine, for instance, that for 
m = mdm and A = X W i m p the relic density constraint is satisfied in the WIMP 
regime, and let us now analyze how the relic density changes as we go from 
A = Xwimp to A = 0. Initially, the relic density will increase, for a smaller A 
implies a smaller annihilation rate and an earlier freeze-out. That is the well- 
known behavior of WIMP particles. It is clear, however, that such behavior 
cannot continue for long because the relic density should go to zero rather than 
to infinity as A — > 0. In fact, if A = the dark matter particle does not interact 
with other particles and cannot be thermally produced in the early Universal. 
Hence, at some point the relic density starts to decrease as A becomes smaller 
and it will continue to do so until it vanishes for A = 0. Since the dark matter 
relic density is a continuous function of A, a smaller value of A, \fi mp , must 
exist such that the relic density is also consistent with the observed dark matter 
density. This alternative solution to the dark matter constraint in the singlet 
model had not been considered before. 

In the FIMP regime, dark matter particles do not reach equilibrium in the 
early Universe and are never abundant enough to annihilate among themselves. 
As a result, the dark matter abundance, Y = n/s, is not determined by the 

2 The corresponding annihilation processes are obtained by exchanging initial and final 
states. 

3 Throughout this paper we do not consider non-thermal production mechanisms. 
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same processes as in the WIMP framework. If dark matter particles are pair 
produced, as it happens in the singlet model, Y satisfies instead the Boltzmann 
equation 

% = ^W M ^)Y eq (T? (2) 

with the boundary condition Y(T » rridm) = 0. In the above equation M p 
is the Planck mass and (av) is the usual thermally averaged production (or 
annihilation) cross section. It is through this quantity that the dependence on 
the particle physics model enters into the evolution of the dark matter abun- 
dance. In the singlet scalar model we have that (crv) oc A 2 . The main difference 
between the WIMP and the FIMP regimes of the singlet model is the value of 
A -or (av)-, which, being much smaller for FIMPs, prevents them from ever 
reaching thermal equilibrium: Y <C Y eq . For that reason, the annihilating term 
proportional to Y 2 does not contribute to the right-hand side of equation © , 
and it cannot be assumed, as is the case for WIMPs, that the dark matter 
particle was in equilibrium, Y(T) = Y eq (T), for T ~ m^,. 

Using equation fl2J) we can anticipate the generic behavior of Y(T) in a large 
class of models including the singlet scalar. Since the right-hand side is greater 
or equal to zero, the abundance either increases or remains constant but never 
decreases, in agreement with the expectation that dark matter annihilations, 
which would reduce the abundance, play no role in this regime. At high tem- 
peratures, Y eq is constant and since all particles are relativistic, we typically 
have that (crv) oc 1/T 2 . In that region then, Y(T) oc 1/T. At low temperatures, 
T « rridm i the particles in the thermal plasma no longer have enough energy to 
produce dark matter particles in their scatterings, (av) — > 0, so the right-hand 
side goes to zero and Y(T) remains constant. Thus, the dark matter abundance 
initially increases as the Universe cools down but at a certain point it reaches 
the freeze in temperature, below which the abundance no longer changes. In 
the next section, we will numerically confirm this simple behavior of the dark 
matter abundance for low and high temperatures. 

Integrating equation from Tj 3> rridm to Tf <C rridm leads to the current 
abundance of dark matter, Y(To). From it, the dark matter relic density can 
be calculated in the usual way: 

flh 2 = 2.742 x 10 8 ^r(T ). (3) 

Notice that the relic density of dark matter in this regime is also determined 
by thermal processes. Besides the given particle physics model, both regimes 
(WIMP and FIMP) assume only the Standard Cosmological Model -e.g. a 
radiation dominated Universe during dark matter production. Regarding dark 
matter, henceforth, they are equally predictive and should be treated on the 
same footing. In the singlet scalar model, there is no objective reason to prefer 
the WIMP solution to the dark matter constraint over the FIMP one. 

Finally, we would like to mention that in spite of the similarity between the 
equations that determine the evolution of the dark matter abundance in the 
WIMP and FIMP regimes, it is usually not possible to compute the FIMP relic 
density using the numerical codes written for WIMPs. MicrOMEGAsjTl], in 
particular, assumes the WIMP paradigm (e.g. Y(T ~ m) = Y eq ), giving wrong 
results if applied directly to the FIMP framework. In our numerical analysis, we 
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Figure 1: The dark matter abundance, Y , as a function of the temperature for different 
values of ms- In this figure mu was set to 120 GeV and A to 10 _ . 



have used micrOMEGAs but only to calculate (av) and have solved equation 
for ourselves. In the next section we present the main results of this paper. 



4 Results 

In this section, we compute the evolution of the dark matter abundance and 
the predicted relic density as a function of the parameters of the singlet model: 
A, ms and nih- We also find the new viable region corresponding to the FIMP 
solution to the dark matter constraint. 

The evolution of the dark matter abundance as a function of the temperature 
is shown in figure[T]for different singlet masses between 1 GeV and 1 TeV. In it, 
A and were fixed respectively to 10~ n and 120 GeV. The general behavior of 

Y is clear: it increases as the Universe cools down until a certain point where it 
freezes in, remaining constant afterwards. As illustrated by the figure, the freeze 
in point depends on the dark matter mass. It is also clear that the rate at which 

Y increases depends on both ms and the temperature. At high temperature the 
singlet behaves as a relativistic particle so Y is pretty much independent of ms, 
as observed in the figure for T > 1 TeV. For ms = 100 GeV and ms = 1 TeV, 
the dark matter abundance freezes in at T ~ ms, when the particles in the 
plasma no longer have enough energy to produce the dark matter particles. 
For ms = 10 GeV and ms — 1 GeV the effect of the higgs resonance, which 
increases (av) significantly, becomes important. That is why Y increases steeply 
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Figure 2: The dark matter abundance, Y , as a function of the temperature for different 
values of A and of ms- From top to bottom on each panel, the lines correspond to 
A = 10 -10 (dash-dotted), A = 10~ n (solid), and A = 10 -12 (dashed). The higgs mass 
was set to 120 GeV. 



for T ~ 100 GeV before freezing in at about 30 GeV. The reason the abundance 
does not keep increasing until T ~ mj for ms = 1 GeV is that, as stated before, 
the production of singlet dark matter is dominated by the W + W~ final state, 
which ceases to be abundant in the thermal plasma well before T ~ 1 GeV. For 
the same reason, the abundances are exactly the same for ms = 1 GeV and 
m s = 10 GeV. 

Let us now see how these results are modified when we vary A and m^. 
Figure [5] shows the dark matter abundance as a function of the temperature for 
four different values of ms (1, 10, 100, 1000 GeV, one in each panel) and three 
values of A: 10~ 10 (upper, dashed-dotted line), 10~ n (middle, solid line) and 
10~ 12 (lower, dashed line). The new feature observed in this figure is that the 
dark matter abundance is, as expected, exactly proportional to A 2 . This is a 
distinctive feature of the FIMP regime: the dark matter yield increases with the 
interaction strength between the dark matter particle and the thermal plasma. 
We will see later that, within this regime of the model, the correct relic density 
of dark matter can be obtained for A between 10~ n and 10 -12 . 

The dependence of the dark matter abundance with the higgs mass is il- 
lustrated in figure [3J It displays Y as a function of the temperature for four 
different values of ms (1, 10, 100, 1000 GeV, one in each panel) and three values 
of the higgs mass: 120 GeV (upper, dashed-dotted line), 150 GeV (middle, solid 
line) and 180 GeV (lower, dashed line). We see that for heavy singlets, which 
freeze in at T > m/i, the abundance does not depend on the higgs mass at 
all. For lighter particles, ms = 1,10 GeV, the effect of different higgs masses 
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Figure 3: The dark matter abundance, Y , as a function of the temperature for different 
values of mh and of ms ■ From top to bottom in each panel, the lines correspond to 
rrih = 120 GeV (dash- dotted) , mh = 150 GeV (solid), and m,h = 180 GeV (dashed). 
In this figure A was chosen to be 10~ n . 



become manifest at low temperatures T < 100 GeV. The smaller the higgs 
mass, the larger the freeze in temperature and the final dark matter abundance, 
as observed in the figure. In any case, since the higgs mass cannot vary over 
a wide range, the corresponding dark matter abundance does not change much 
either with m^. Its total effect amounts to a variation of less than one order of 
magnitude in Y. 

In the FIMP regime, as the previous figures illustrate, the dark matter par- 
ticles are slowly produced as the Universe cools down until the freeze in temper- 
ature is reached. In the singlet model, the freeze in temperature is given by the 
singlet mass for ms > 100 GeV and by T ~ 30-50 GeV for lighter singlets. Close 
to the higgs resonance the production accelerates but the abundance remains 
small (Y <C 1) throughout the entire history of the Universe. In the FIMP 
regime, the dark matter particles never reach thermal equilibrium or annihilate 
among themselves. 

Now that we have analyzed the evolution of the dark matter abundance in 
the FIMP regime, we turn our attention to the relic density, the quantity that 
is actually constrained by observational data. The present relic density of dark 
matter is proportional to ms and to the asymptotic value of Y, see equation 
Figure [H displays the predicted relic density as a function of the singlet mass 
for m/j = 120 GeV and different values of the coupling A. From top to bottom 
the lines correspond to A = 3 x 10~ n (dashed line), A = 10~ n (solid line), 
and A = 3 x 10~ 12 (dashed-dotted line). For comparison, the region compatible 
with the Wmap determination of the dark matter density [T5] is also shown as a 
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Figure 4: The relic density of dark matter, Q.h 2 , as a function of ms for different 
values of X. The higgs mass, mn, was fixed to 120 GeV in this figure. The Wmap 
compatible region is also shown as a horizontal band. 



horizontal band. As anticipated, the relic density is proportional to A 2 . Notice 
that for A = 3 x 10~ n (the dashed line) the relic density is, independently of 
ms, too high to be compatible with the Wmap range . For the other two values 
of A shown, we do find regions consistent with the data. Typically, in the singlet 
model there are two different masses that are consistent with the dark matter 
constraint for a given value of A. It is clearly observed in the figure that the 
relic density has a sizable jump at the higgs resonance, ms ~ 60 GeV, or more 
generally at ms ~ m^/2. If ms is slightly below the resonance the relic density 
is significantly larger than if it is slightly above the resonance. Moreover, while 
above the resonance the relic density is essentially independent of the singlet 
mass, it becomes proportional to ms below the higgs resonance. That is exactly 
the region where we have found the abundance, Y, to be independent of ms 
-see previous figures. 

The variation of the relic density with the higgs mass is shown in figure [5] 
We can see that the region where the behavior of the relic density changes is 
indeed determined by the higgs resonance condition: ms ~ rrih/2. For high 
masses, mj > m^, the relic density is seen to be independent of m^. For lower 
masses, ms < rrih/2, the relic density turns out to be about a factor 3 larger 
for m/i = 120 GeV than for rrih — 180 GeV. Hence, in that region, the value 
of the mass that is compatible with the relic density might change by the same 
factor as m/i varies. 

Usually we are interested not in computing the relic density by itself but 
rather in obtaining the viable parameter space of the model, the regions where 
the predicted relic density is consistent with the observed dark matter density. 
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Figure 5: The relic density of dark matter, Q.h 2 , as a function of ms for different 
values of rrih- In this figure the coupling A was set to 10 -11 . The Wmap compatible 
region is also shown as a horizontal band. 



Figure H2 shows, in the plane (ms,A), the viable region of the singlet model 
for the FIMP regime (see e.g. [12] for the corresponding regions in the WIMP 
regime). Each band shows the region compatible with the dark matter con- 
straint (at 2a) for a given higgs mass. Notice that for heavy dark matter, 
ms > 100 GeV, A should be about 1.2 x 10~ n independently of the singlet 
mass or the higgs mass. For dark matter masses between 0.1 GeV and 50 GeV 
the required A decreases from about 2 x 10~ n to about 10 -12 , reaching this 
minimum value for ms ~ 40 GeV and mt = 120 GeV. Around the higgs reso- 
nance, the value of A that satisfies the dark matter constraint varies by about 
one order of magnitude in a narrow range of mg. Notice if ms < 100 GeV for 
any given value of A there usually exists two solutions for ms- If, for instance, 
m/j is found to be 150 GeV, for A = 3 x 10~ 12 the dark matter mass could be 
either about 8 GeV or about 75 GeV. In any case, the crucial point is that such 
solutions do exist. 

This last figure clearly demonstrates our main result, that in the singlet 
model the dark matter constraint can also be satisfied in a previously unexplored 
region of the parameter space. In this new viable region we have unveiled, 
the interactions between the singlet and the Standard Model fields are very 
suppressed and the singlet behaves as a Feebly Interacting Massive Particle 
(FIMP) rather than as a WIMP. We would like to emphasize that this FIMP 
regime of the singlet model is as natural and predictive as the WIMP regime. 
There is no need to include new particles, to introduce non-thermal production 
mechanisms, or to modify the cosmological model in order to satisfy the dark 
matter constraint. From the dark matter point of view, therefore, both regimes 
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Figure 6: The viable parameter space of the singlet model in the FIMP regime. Inside 
the bands the dark matter constraint is satisfied for the given value of the higgs mass. 



should be treated on equal footing. 



5 Discussion 

A generic property of the new viable region we have found is the absence of any 
direct or indirect detection signature of dark matter, a direct consequence of 
the feeble interactions of the dark matter particles. This lack of experimental 
signatures, which may currently seem as a drawback, could become a desirable 
feature if in the next years dark matter is not detected. Such an event would 
force us to abandon the WIMP paradigm and to replace it with an alternative 
framework. Among those that have been considered so far, FIMP dark matter- 
looks particularly attractive due to its simplicity and predictivity. 

Even though the FIMP regime of the singlet model cannot be directly probed 
in dark matter experiments, it can certainly be falsified. As soon as a credible 
signal of dark matter is observed in direct or indirect detection experiments, the 
FIMP regime of the singlet scalar model can be excluded as the explanation of 
dark matter. 

Recently, it was found that the latest constraints on dark matter direct 
detection by the XENON100 experiment [TB] exclude a significant fraction of 
the parameter space of the singlet model -see figure 3 in [17 . In particular, 
tos < 80 GeV seems to be strongly disfavored. We would like to stress, however, 
that such a bound applies only to the WIMP regime of the singlet model and that 
no comparable bound exists in the FIMP regime. In it, the region ms < 80 GeV, 
as well as the entire GeV to TeV mass range, is perfectly compatible with the 
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dark matter constraint and with the absence of a dark matter signal in direct 
and indirect dark matter experiments. 

Finally, we would like to mention that the analysis we have presented for 
the singlet scalar model might be extended to other models of dark matter. As 
it is clear from our previous discussion, the main requirement on such models 
is that the interactions of the dark matter particle depend on a free parameter. 
If that is the case, it is likely that in such models also exist a new dark matter 
compatible region where the relic density constraint is satisfied via freeze in of 
Feebly Interacting Massive Particles. 

6 Conclusions 

The singlet scalar model is one of the simplest extensions of the Standard Model 
that can explain the dark matter. It contains just one additional field, which 
is neutral under the SM gauge group and odd under a new Z 2 symmetry, and 
it introduces only two free parameters: the dark matter mass and the coupling 
between the dark matter and the higgs, A. We have demonstrated that in this 
model there exists a previously unexplored region of the parameter space where 
the dark matter constraint can naturally be satisfied. In it, the dark matter be- 
haves as a Feebly Interacting Massive Particle, FIMP, and its relic density is the 
result of a thermal freeze in in the early Universe. We have examined in detail 
the FIMP regime of the singlet scalar model. After reviewing the calculation 
of the relic density, the evolution of the dark matter abundance with the tem- 
perature was analyzed as a function of the parameters of the model. Then, we 
computed the dark matter relic density and obtained the new viable parameter 
space of this model. It features couplings of order 10 -11 -10 -12 for dark matter 
masses in the GeV to TeV range. Finally, we discussed the implications of this 
new viable region for the correct interpretation of recent experimental bounds. 
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